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ABSTRACT 



We have analyzed all the observations of BP Tauri taken by the International 
Ultraviolet Explorer in the low resolution (AA ~ 6A), long wavelength (from 
A = 1850A to A = 3350A) range. This dataset contains 61 spectra. We observe 
variability in the ultraviolet continuum of Am cont . ~ 1 magnitude and variability 
in the Mg II line flux of AmM g n ~ 0.8 magnitudes. Moreover, these spectra do 
not show any correlation between the continuum flux and the Mg II line flux, 
thus resolving a standing controversy in the literature concerning the origin of 
the Mg II line flux. There is no correlation between the color temperature of the 
UV continuum and the average value of its flux. Using models of the accretion 
process developed by Calvet & Gullbring (1998), we obtain energy fluxes, accre- 
tion spot sizes, and accretion rates from the IUE observations of BP Tauri. We 
find average energy fluxes of 5.0 10 11 ergs cm -2 s _1 , average spot sizes of 4.4 10~ 3 
times the stellar surface, and average accretion rates of 1.6 lO~ 8 M /yr. Our 
analysis shows that the particle energy flux and the UV flux in the stellar surface 
are proportional to each other. Most strikingly, we observe a correlation between 
accretion rate and spot size, with the spot size increasing as the square of the 
accretion rate. Based on the results of a simulation, we conclude that geometrical 
effects (i.e. the varying inclination of the spot with respect to the observer) are 
not enough to account for this effect. Current models of the accretion process fail 
to reproduce such an effect, suggesting the need of using more realistic descrip- 
tions of the stellar field when treating magnetospheric accretion. There may also 
be an unmodelled efficiency factor that determines how matter is loaded into 
the field lines. Non-dipole fields, geometry, oblique shocks and the possibility 
of "limb brightening" should be taken into account when creating models and 
explaining the results of observations of T- Tauri stars. 
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Subject headings: stars: pre-main-sequence — stars: individual (BP Tauri) - 
space vehicles — ultraviolet: stars 

1. Introduction 

The International Ultraviolet Explorer (IUE) was operational for almost 20 years 
and took over 107000 images of a wide array of objects. In particular, it observed about 
130 T- Tauri stars (TTSs) in the low resolution (AA ~ 6A), long wavelength range (from 
A = 1850A to A = 3350A). This dataset constitutes the most complete set of observations 
of TTSs in the Balmer range. 

Now that the IUE Final Archive has been completed and all the data is in the public 
domain 1 we have undertaken the project of exploring all this dataset. In this paper we will 
examine the observations of BP Tauri (also known as HBC 32 and HD 281934). A complete 
discussion of the full dataset will appear elsewhere. 

We have chosen this star in particular because of the large number of observations by 
IUE: there are 61 low resolution spectra. They will allow us to study correlations between 
different physical aspects of the T-Tauri phenomenon over a long period of time. 

BP Tauri is a single (Ghez, Neugebauer & Matthews 1993), Classical T-Tauri star 
(CTTS) located in L1495 (Lynds 1962) in the Taurus-Auriga complex of dark clouds at a 
distance of 140 pc (Elias 1978). The star has a constant radial velocity close to that of its 
molecular surroundings and a projected equatorial velocity of i>sim < 10 km/s (Hartmann 
et al. 1986; Basri & Batalha 1990). The spectral type of BP Tauri corresponds to a K5-K7 
dwarf with a strong absorption in the lithium resonance line (Duncan 1991; Hartigan et al. 

x http: //archive, stsci.edu/iue/ 
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1989). The mass is M = 0.5 - 0.8M© (Gulbring 1994; Gullbring et al. 1998). 

Photometric observations (Vrba et al. 1986; Bouvier, Bertout & Bouchet 1988; Simon, 
Vrba & Herbst 1990 - SVH90 - and the observations summarized by Rydgren et al. 1984) 
have shown that BP Tauri has a large infrared emission and is variable in all optical bands, 
with a range in B of up to 3 magnitudes (Herbig & Bell 1988). Rapid fluctuations on 
time-scales of minutes have been observed by Schneeberger, Worden & Africano (1979). 

These fluctuations make it difficult to extract a possible photometric period. No period 
at all (Bouvier et al. 1988, Simon et al. 1990) and periods of 6.1 days (Simon et al. 1990), 
7.6 days (Vrba et al. 1986), 7.7 days (Osterloh, Thommes & Kania 1996), 8.3 days (Richter 
et al. 1992), and ~ 10 days (Gomez de Castro & Franqueira 1997 - GdCF97) have been 
proposed. The fact that BP Tauri grows redder as it becomes fainter has been used to 
explain this variability in the period in terms of a hot spot on the surface of the star (see for 
example SVH90). In a small fraction of CTTSs, spots cooler than the photosphere have also 
been observed, but no evidence of cool spots exists for BP Tauri (Gullbring et al. 1996). 
Variations in the measured period of BP Tauri might correspond to the hot spot (or spots) 
moving to different latitudes in a differentially rotating star. An inclination of i ~ 30° — 50° 
for BP Tauri has been derived on the basis of the measured v sin i and the longest measured 
periods (SVH90; Gullbring 1994). The variability in the photometric period can also be 
explained by considering a beat frequency between the Keplerian frequency of clumps of 
material circling the star and the rotation rate of the stellar surface (Bouvier et al. 1999; 
Smith, Bonnell & Lewis 1995; Smith, Lewis & Bonnell 1995). 

Measurements of spot sizes in this and other TTSs give values that vary between a 
fraction of a percent to almost 40% the visible stellar surface (Vrba et al. 1986; Vrba et al. 
1993; Herbst, Herbst & Grossman 1994; Gullbring 1994; Bouvier et al. 1995; Fernandez & 
Eiroa 1996). 
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As has been observed in many CTTSs, the spectrum of BP Tauri appears veiled by 
variable continuum excess emission. This is the origin of the photometric variability. The 
optical veiling has been observed to be between 30 and 75% the value of the photospheric 
flux (Basri & Batalha 1990; Valenti, Basri & Johns 1993; Hartigan, Edwards & Ghandour 
1995; Gullbring et al. 1998). 

Some spectral lines (notably, the Balmer series and the Mg II line) show evidence of 
blueshifted and redshifted absorption (see for example, Johns & Basri 1995ab; GdCF97). 
It has also been noted that line emission and the continuum level are sometimes correlated 
(GdCF97; SVH90). We will explore one of these correlations in Section 3. 

These observations can be framed in terms of the current paradigm for CTTS. As 
we understand it today, a CTTS consists of a late-type pre-main-sequence star with a 
magnetic field, surrounded by an accreting, dusty disk. At the distance from the star at 
which the torque due to the magnetic field equals the viscous torque in the disk, the disk 
is truncated. Some material from the disk will then be captured by the magnetic field (see 
for example Konigl 1991) and some will escape as a wind (Shu et al. 1994). Collisional 
coupling between ions and neutrals is such that even the neutral gas component follows the 
magnetic field lines (Martin, 1996). At the stellar surface, the material will have supersonic 
velocity. The simplest models assume that this velocity should be the free-fall velocity 
(vff ~ 300 km s -1 , Calvet & Gullbring 1998) but more complicated models slow down the 
material to a fraction ot this velocity (Ostriker & Shu 1995). Nevertheless, a shock will 
occur some distance above the star. The shocked gas will heat up the surface of the star 
and a hot spot will be seen. Emission from this hot spot is thought to be responsible for 
the visible and ultraviolet excess observed in CTTS spectra (Valenti et al. 1993; Calvet 
& Gullbring 1998; Lamzin 1998). Furthermore, emission from the whole accretion flow is 
believed to be responsible for other features of the TTS spectra, in particular, the strong 
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emission in the Balmer lines (Hartmann, Hewett & Calvet 1994; Johns-Krull & Basri 1997) 
and some atomic lines (Batalha et al. 1996; Beristain, Edwards & Kwan 1998). In this 
model redshifted absorption features would correspond to material falling into the accretion 
spot. Blueshifted absorption features would correspond to absorption by an optically thin 
wind. It is unclear at this point what fraction of the excess emission (if any) is due to the 
wind. 

In the context of the magnetospheric accretion model, the ultraviolet is a particularly 
interesting range. Emission from the underlying star is low compared to that of the veiling 
continuum. This means that, even though is not possible to measure the veiling directly (as 
it is possible in the visible range), conclusions derived from the ultraviolet analysis are not 
very sensitive to errors in the spectral type. 

In this paper we will explore the IUE dataset for BP Tauri. Relatively little work has 
been done on the individual IUE observations and their relation to the accretion paradigm. 
Here we assume the excess ultraviolet continuum emission to be caused by accretion. We 
start with a discussion of the raw data. A more complete discussion (including many more 
stars) will appear soon, in a paper dealing with the full low-resolution IUE dataset. In 
Section 3, we will use our dataset to study the correlation between the UV continuum and 
the flux in the Mg II blend. To understand the IUE data to the fullest we need a physical 
model. Therefore in Section 4 we will use the Calvet & Gullbring (1998) (CG98) description 
of the accretion shock to interpret our spectra in terms of accretion rates. Section 5 and 6 
contain the analysis and the conclusions. 
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2. Observations 

We searched the IUE Merged Log for observations of TTS taken with the long 
wavelength cameras (LWP and LWR, from A = 1850A to A = 3350A) and obtained spectra 
for a total of 131 different TTS. For BP Tauri we found 61 low resolution (AA ~ 6A) 
spectra. A summary of the observations is shown in Table 1. The oldest spectra comes 
from April 1979 and the most recent from January 1992. There are clusters of observations 
taken in October 1986 (SVH90) and January 1992 (GdCF97). 

All the spectra used in this paper have been processed by the New Spectra Image 
Processing System (NEWSIPS). For an explanation of this system and the differences with 
IUESIPS, see Nichols & Linksky (1996). For each observation the spectrum has been 
sky-subtracted by using the background in the camera. NEWSIPS spectra come with an 
estimate of the error at each wavelength. This error is derived empirically for each camera 
by measuring the scatter in the flux around the mean in the background regions of several 
hundred images (Garhart et al. 1997). As such, it is an estimate of the instrumental, 
statistical and background noise for each spectrum. NEWSIPS also provides a "quality" 
flag at each wavelength of each spectrum. This flag indicates certain special events, such as 
a saturated pixel, cosmic ray or the presence of a calibration mark. 

Figure la. shows the mean and the normalized variance (Johns & Basri 1995b) of all 
spectra. The feature at A ~ 3050A is a reseaux mark. The central line is the blending of 
the Mg II h & k lines. Figure la shows that the edges of the spectra vary the most, whereas 
the center is very constant. The noise of the spectra also increases towards the edges (due 
to the reduced sensitivity of the Vidicon detectors used by IUE) but not enough to account 
for the value of the variance. 

Figure lb shows two extreme manifestations of the spectra (not corrected for 
reddening). The lower one (lwpl8976, solid line) shows strong emission of C II] and Fe II 
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multiplets. For a more exhaustive description of the spectral features present, see Imhoff & 
Appenzeller (1987) and Brown, de M. Ferraz & Jordan (1984). In the upper one (lwp09256, 
dash-dot line), those emissions are very difficult to distinguish. These spectral changes are 
typical of CTTS and are probably due to inhomogeneities in the accretion flow and/or 
changes in the orientation of the star as it rotates. The temperature of the upper spectrum 
is almost a thousand Kelvin higher than that of the lower spectrum (Table 1). 

We have de-reddened all the spectra using Ay = 0.51 magnitudes (Gullbring et al. 
1998). This reddening determination is based on the shape of the optical spectrum and not 
only on broad band photometry. It is therefore more reliable than previous values of Ay. 

In addition to its science objects, IUE also observed a suite of standard stars. We 
have compared spectral characteristics of BP Tauri with two of these standards (Table 1). 
The color temperature of the BP Tauri spectra was evaluated using the continuum fitting 
procedures described in Sections 4.2 and 5.3. HBC 399 (also known as V 827 Tau) was the 
only WTTS observed by IUE with enough signal-to-noise to measure its color temperature. 
It is classified as K7, MO and has an Ha: equivalent width of — 1.8A (Herbig & Bell 1988). 
We take the reddening to be Ay = 0.28 magnitudes (Kenyon & Hartmann 1995). IUE 
also observed main-sequence stars including 61 Cyg B, with spectral type K7V. According 
to the SIMBAD database its parallax is ~ 0.28 arc-sec, which means that the reddening 
is negligible. As Table 1 shows, the continuum level and the flux in the Mg II line in BP 
Tauri are one and three orders of magnitude larger than those in HBC 399 and 61 Cyg B 
respectively. This observation confirms the fact that the intrinsic photospheric continuum 
of the star is invisible in the Balmer range. As shown in Section 4 this ultraviolet excess 
can be understood as being due to the accretion process. 
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3. The continuum and the Mg II blend 

As mentioned in the introduction, relationships between line emission and continuum 
level have been published by a number of authors. The study of these relationships is 
interesting because they illuminate what fraction of the observed line emission is due to 
the accretion process, as opposed to the wind or the stellar photosphere. For BP Tauri in 
the ultraviolet, two main articles reaching opposite conclusions have been published. These 
articles use subsets of the entire IUE dataset for BP Tauri. 

From six spectra taken in 1986, SVH90 reached the conclusion that variations in the 
flux of Mg II blend can be correlated with continuum variations. They also looked at O I 
(1300A), Si II (1810A) and C II (1335A) and reached the same conclusion. This would 
seem to indicate that the bulk of the emission in those lines comes from a process directly 
related to the formation of the accretion spot. From a set of twenty- two spectra taken in 
1992, GdCF97 reached the conclusion that the Mg II blend, Si II and C IV (1550A) light 
curves are dominated by short timescale fluctuations, whereas O I and He II (1640A) are 
well correlated with variations in the UV continuum. Mg II, Si II and C IV are collisionaly 
excited lines, but O I and He II can be excited by recombination process. 

The IUE data utilized in this article includes all the data used by SVH90 and GdCF97 
and much more and therefore it can resolve the issue of whether or not the Mg II line 
flux correlates with the UV continuum. Following GdCF97, we have measured the UV 
continuum variations in a window of 50A width centered at 2900A. This region of the 
spectra stands between two groups of Fe II lines (Brown et al, 1984) and therefore it 
samples closely the real continuum of the star. 

Disregarding those spectra for which the continuum or the Mg II blend are affected by 
"extraordinary" events (defined as having a "quality" flag different from zero), we have 43 
measurements. The dispersion in the continuum flux measured in this way, is Am cont . ~ 1 
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magnitudes. The dispersion in the flux of the line is AmMgii ~ 0.8 magnitudes. Figure 
2 shows that, indeed, no correlation exists between the flux in the blend and the UV 
continuum flux. 

The difficulty of interpreting low- resolution observations of the Mg II line is compounded 
by the fact that high resolution observations, as described by GdCF97 show that the line 
has a strong wind component. If part of the line is created in the accretion flow and part in 
an extended wind, a measurement of the low resolution flux is sampling to different regions 
of the matter flow (see Section 5.3). 

We have also compared the color temperature of the continuum (with the continuum 
as found in Section 4.2) and the amount of continuum flux (Figure 3). The errors in the 
temperature are ~ 10%. As Figure 3 shows, no correlation exists between color temperature 
and flux in the continuum. This implies that there is no correlation between the flux in the 
line and the color temperature. We defer further discussion of this Figure until Section 5.3. 

4. Modeling the ultraviolet continuum 

In order to interpret the spectra, a physical model is needed; we have decided to use 
the work of CG98. A full description of the model can be found there. We describe only 
those aspects of it directly relevant to our data. In this section we will find the accretion 
rate, the energy flux and the spot size indicated by each of our spectra. 

4.1. The Model 

CG98 construct a model in which the observed spectral excess in CTTS is the sum of 
optically thick emission coming from the heated photosphere, and optically thin emission 
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coming from the pre- and postshock regions. The shock is assumed to be parallel to the 
surface of the star. The heated postshock material radiates mainly in X-rays. Half of that 
radiation is reprocessed into optical and UV wavelengths in the photosphere; the other 
half is absorbed in the preshock region and reprocessed in optical and UV wavelengths. A 
fraction of this emission will be re-emitted towards the star and will contribute to heating 
the photosphere, whereas the remainder will travel away from the star. 

The amount of observed radiation at each wavelength depends on two factors. The 
first factor is the strength of the potential depth of the star, M/R, where M is the mass 
and R is the radius of the star. This ratio determines the temperature of the postshock gas 
and therefore, everything else being equal, the hardness of the incident radiation. We have 
assumed M/R = O.25M /R for BP Tauri (Gullbring et al. 1998). The second factor is the 
kinetic energy flux carried by the accretion flow. From mass conservation we have that the 
density of the accretion flow is p = M / (A int _ spot * f //), where M is the accretion rate, v ff is 
the free-fall velocity of the material. A int - S pot is the intrinsic area of the spot, equal to the 
observed area divided by cos#, the inclination of the spot with respect to the line of sight. 
To consider inclination effects explicitly, we write p = Mcos9/(A spot * f//), where now A spot 
is the observed spot area. We then write the kinetic energy flux (or the particle flux) as: 

^-1 i. -^m/ McOS# w M w it! . o. f , , o ^ /-.x 

JT = -pvff 3 = 9.8 10 10 -)( ) (— ergs cm" 2 s" 1 1 

2 H " V 1O" 8 M yr- lJK 0.5M Q JK 2R & J v 0.0r B KJ 

in which f// — 300 km sec -1 , the free-fall velocity from a truncation radius equal 
to Rt = 5R* (Meyer, Calvet & Hillenbrand 1997), / is the apparent fractional size of 
the accretion region or spot, A spot = f4:iiRl and cos 6* is the inclination of the spot with 
respect to the line of sight. The factor of cos# is obtained by assuming that the star has 
only one spot. If the region producing the emission is not one spot but a set of spots or a 
deformed ring (Mahdavi & Kenyon 1998) cos# should be interpreted as the cosine of the 
mean inclination of the region responsible for the emission. In the following analysis we do 
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not find the intrinsic accretion rate but Mcos#, the "projected" accretion rate. Similarly, 
the value that is found for /, the spot size, corresponds to the observed spot size and not 
the intrisic one. Both are related by f — fi nt cos 9. 

CG98 assumes that the truncation radius Rt remains constant. Modeling by a number 
of authors (e.g. Konigl 1991, Ghosh 1995, Ostriker & Shu 1995) suggests that R T a M~ 2 / 7 . 
This result can be obtained in a number of ways, but the simplest is by realizing that it 
is the only dimensionally correct dependence between Rt and M. Furthermore, in most 
models the truncation radius affects the size of the spot in such a way that a large Rt 
implies a small spot (see Section 5.6). If we were to assume a varying truncation radius, 
the factor ( = 1 — R*/ R T would appear multiplying the accretion rate. As the truncation 
radius depends only weakly on the accretion rate, we will assume ( to be constant . 

CG98 divides the emission accretion zone into three regions: Shock/postshock, 
preshock, and heated atmosphere. In the shock and postshock region a number of 
simplifying assumptions are made to solve the mass, momentum, and energy conservation 
equations: the flow is assumed to be plane-parallel and strictly along field lines (i.e., 
one-dimensional) perpendicular to the stellar surface. Energy transport by conduction is 
neglected, as are instabilities on the flow and effects of the magnetic field on the dynamics of 
the gas. For the boundary conditions the shock is assumed to be strong. The temperature 
structure towards the stellar surface is calculated until a minimum is reached. The emission 
from the accretion column is calculated with CLOUDY. 

The preshock region is assumed to be illuminated by the radiation field provided by 
the postshock calculation. The small size of this region compared to the stellar radius is 
used to justify an assumption of constant density through it. 

The stellar atmosphere is assumed to be semi-infinite, plane-parallel, and in radiative 
equilibrium. The emission from it is the sum of the unperturbed stellar photosphere, plus 
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reprocessed incident radiation. For T* = 4000 K, M = O.5M , R = 2R and logJF = 11.5, 
the effective temperature of the continuum at r Ross = 2/3 is then T cff ~ 6000K. The output 
of the model is a grid of the excess spectra each of which is characterized by a different 
value of log T . 

Through CG98's model, the excess Balmer continuum is shown to be due to optically 
thick emission from the heated atmosphere and optically thin emission from the pre- and 
postshock regions. The former dominates at high energy fluxes (logjF > 11) with the total 
spectrum becoming more like blackbody emission. The optically thin emission becomes 
more important as the energy flux decreases. The Paschen continuum is always dominated 
by optically thick photospheric emission. We refer the reader to the original paper for 
further details. 

4.2. Fitting the continuum 

To use CG98 models, the first step is to find the UV continuum. We have only used IUE 
spectra in which there are regions with S/N>2.5. As this region differs for each image, the 
useful wavelength coverage of each spectrum varies slightly, but in general, the useful range 
is reduced from the nominal IUE coverage to A ~ 2400 A to A ~ 3 100 A. The signal-to-noise 
constraint reduces the number of useful spectra to 45. We have developed an automated 
procedure that eliminates from each spectrum the regions with the resonant Mg II blend 
and the Fe II uv 1 blend (2771.3 to 2821.3 A and 2576 to 2636 A, respectively). For the 
remaining data, the bottom of the flux in a AA ~ 10A wide box is found. This defines the 
continuum. The error in the continuum determination is found from the variations in the 
flux within this box. The continuum points found with this procedure are binned until we 
have 20 data points. 
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As Figure lb shows, there are certain regions of the spectrum that sample the apparent 
continuum. Our procedure makes sure that the continuum spline fit goes through certain 
points like 2300A, 2650A, 2900A and 3000A. 

4.3. The Temperature of the Underlying Star 

To obtain the spectrum of the excess emission we must first characterize the spectrum 
of the underlying K7 star. As shown in Table 1 the color temperature of the standard stars 
of the same spectral type as BP Tauri are lower than the canonical 4000 K of a K7 star. 
This is because the UV photosphere is optically thicker than the visible one and therefore 
it is situated close to the temperature minimum of the star. If we assume that BP Tauri 
is intrinsically like HBC 399, the color temperature of its photosphere should be ~ 3200K. 
It has been shown that the ultraviolet color temperature of early G stars is ~ 80% of their 
effective temperature (Haisch & Basri 1985). It is also shown by CG98 that the ratio of 
Balmer to Paschen continua for a T cff = 4000K star is better reproduced by assuming that 
the underlying star has a temperature of 3580 K in the Paschen continuum. They attribute 
this to the possibility that accretion is occurring into a cool spot. What the temperature of 
this spot is in the Balmer range cannot be estimated without further modeling. 

Therefore, for the purposes of calculating the excess emission, we assume that the 
underlying star has a temperature of 3580K. This estimate is consistent with the color 
temperature of the standards and allows us to use CG98 models. We will examine the 
consequences of this assumption in Section 5.4. 
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4.4. Finding the Energy Flux and the Spot Size 

Once the underlying star has been subtracted, the resulting observed excess emission 
is compared with model spectra given by CG98. To perform this comparison, we obtain 
the mean value and the slope (assuming a straight line) of the logarithm of the observed 
spectrum. These two parameters have zero covariance and so any correlation between them 
is not due to the fitting procedure. 

In CG98 models, the slope of the predicted excess emission for a fixed value of M/R 
depends only on log T . By fitting the continuum we are able to associate a value of log T 
with each observed spectrum. When our values for \ogT extend beyond their published 
grid, we have linearly extrapolated their models. Once the value of \ogT is found, we 
find the relative size of the spot / by comparing the mean value of the theoretical excess 
spectrum to that of the observed excess spectrum. 

In this way we obtain the values for \ogT and / quoted in table 2. Errors in / are 
<jf/f ~ 10% and are a combination of the errors in the value of the flux at each wavelength, 
and errors in the determination of the continuum. Errors in log J 7 are a\ og jr /log J 7 ~ 1%. 
With Equation 1 we obtain Mcos9, the accretion rate uncorrected for projection effects. 
With this model it is impossible to obtain the "unprojected" accretion rate M, unless cos 9 
is known: T is independent of 9 but the observed spot area / is proportional to cos 9 and so 
their product is proportional to M cos9. The errors in M cos 9 are cos9 / M cos9 ~ 30%. 
All errors include an error in Ay of 0.05 magnitudes. 

5. Analysis 

The average size of / is / = 4.4 10~ 3 with a dispersion of a — 5 10~ 3 . The average 
of logF is log J 7 = 11.7, and the dispersion is a = 0.2 (see Figure 4a and 4b). It should 
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be noted that the value of / is a lower limit to the real spot size, given that what one 
observes is the projected spot area. Also, because of the period changes mentioned in the 
introduction, we may not be observing the same spot (or set of spots) all the time, but their 
size and distribution over time. Figure 4c shows the distribution of Mcos9. The average 
is Mcos9 = 1.6 1CT 8 M yr -1 , and the dispersion is a ~ 0.8 lO~ 8 M Q /yr. The dispersion 
of / is larger than the observational errors, but the dispersions in log T and M cos 9 are 
consistent with the errors. 

As mentioned in the introduction, a wide range of spot sizes have been identified by 
other authors. We find spot sizes ranging from 0.06 to 3% the size of the visible star disk. 
The average accretion rate we observe agrees with the most recent determinations for BP 
Tauri, based on single optical spectra: M ~ 2 — 3 1O~ 8 M yr _1 (Gullbring et al. 1998, 
CG98). 

The logarithm of the energy flux, logjF, does not depend on the inclination angle 
9. Therefore, the variations observed in Figure 4b are due to variations in M and/or the 
intrinsic size of the accretion spot. The dispersion in log J 7 implies a fractional variation 
of ~ 50% in M/fint, where f int is the intrinsic (i.e. deprojected, f int = f / cos 9) spot size. 
This dispersion is consistent with the errors in T. Variations in M cannot be decoupled 
from those in f int . What can be said is that changes in the accretion rate M or f int are less 
than 50% during the time spanned by the observations. This upper limit puts constraints 
on the inhomogeneities of the disk (that may contribute to the variations in the accretion 
rate) and in the geometry of the stellar field. 
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5.1. Variable Extinction 

Recently, Bouvier et al. (1999) reported on photometric observations of the CTTS 
AA Tauri that can be interpreted as periodic occultations of the stellar photosphere by a 
warp in the accretion disk. One of the effects of this warp is to act as a source of variable 
reddening. This mechanism might also play a role in the IUE observations of BP Tauri. To 
find out if that is the case, we performed the following analysis. The procedure of finding 
which model corresponds with an observed spectrum has three free parameters: reddening, 
/ and log J 7 . In this paper we assume a fixed reddening and find the values of / and logjF 
that match the observed spectrum to the models. We can also assume a fixed spot size 
and find the values of log T and reddening that match the observed spectra. We did so, 
assuming / = 4 1CT 3 . For the cases in which this procedure succeeded in finding a logT and 
A v that fit the data, we obtained a large range of values for A v , from 0.4 to 10 magnitudes, 
but most of them bigger than 1. The largest value that we found in the literature for a 
reddening determination of BP Tauri was 0.85 (Valenti et al. 1993). Therefore, most of the 
reddenings determined in this way do not fit the optical data. Photometric observations in 
various bands simultaneously have also failed to detect color changes in BP Tauri analogous 
to those reported by Bouvier et al. (1999) for AA Tauri. We therefore conclude that the 
changes in log J 7 and / are not due to changes in reddening. 



5.2. The accretion luminosity versus the continuum flux 

Figure 5 shows a tight correlation between Mcos9 and the continuum flux, which 
implies that the observed (projected) accretion rate determines the brightness of the UV 
continuum. We can convert the observed UV flux to luminosity and the model accretion 



rate to accretion luminosity. We obtain 

log(L acc /L ) = (1.1 ± 0.1) * log(L 7L ) 



- (0.5 ±0.05) 



(2) 



A similar relationship has been found by GC98 and Gullbring et al. 1997, who observe 
a very clear correlation between the accretion luminosity and the brightness in the U-band, 
based on individual observations of various CTTSs in the optical. The empirical relation 
they find is 

log(L acc /L ) = 1.09 * \og(Lu/L e ) + 0.98 (3) 

in which L acc is the accretion luminosity (which is proportional to the accretion 
rate) and Ljj is the unreddened luminosity of the star in the U-band once the intrinsic 
(photospheric) luminosity has been subtracted. 

An independent way to verify the accretion rates obtained here is to use Equation 3 as to 
predict accretion luminosity and then compare it with the accretion luminosities from Table 
4. Given the rapid photometric variations of CTTSs, one needs to use concurrent L acc and 
Ljj to test this relationship. We searched the photometric database maintained by William 
Herbst 2 described in Herbst et al. (1994) and found two measurements of the U-magnitude 
close enough to the IUE measurements: RJD=5615.51 (IUE at 5615.85) and RJD=5620.99 
(IUE at 5620.97), which imply \og(L v / L & ) = -1.88 and \og(Lu/L e ) = -2.05 respectively. 
These determinations assume LspTauri = O.925L (Gullbring et al. 1998). Using Equation 
3, the predicted values are log(L acc /L ) = —1.072 and log(L acc /L ) = —1.25 respectively. 
We observe log(L acc /L ) = —1.101 and log(L acc /L ) = —1.303, by interpolating Figure 5. 
This is consistent with Equation 3 within the errors of our data. Admittedly, two points 
out of 45 is a very limited sample, but the fact that they obey Equation 3 lends support to 
our determinations of the accretion rate. 



2 http://www. astro. wesleyan.edu/ bill/research/ttauri.html#ov2 
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5.3. The Line Flux and the Continuum Temperature 

Figure 3 shows that no correlation exists between the continuum temperature and the 
flux in the continuum. Within the context of CG98 model it is possible to understand 
this result. The Balmer continuum is the sum of optically thick emission from the heated 
photosphere and optically thin emission from the pre-shock region. Only when the value of 
T is very high does the observed color temperature approach the temperature of the heated 
photosphere. The flux in the continuum at 2900A depends, again in the context of CG98's 
model, not only on the spot size, but also on the value of T . Two spectra with the same 
accretion spot size but different T will have different values of continuum flux. Therefore, 
the color temperature is a measure of two different things, the spot size and the energy flux. 

We do not observe any correlation of the Mg II line flux with log T, the accretion rate 
Mcos#, or the spot size /. 

It has been argued by Hartmann (1990) that strong emission lines (like Ha and Mg II) 
in the spectra of CTTS are formed in an extended envelope. His argument is based on the 
similar widths of the strong lines. Table 1 shows that the Mg II line is very bright in BP 
Tauri, compared with that of non-CTTS of the same spectral type. There have also been 
unsuccessful attempts to model the Mg II line blend as purely chromospheric (e.g. Calvet, 
Basri & Kuhi 1984). All these arguments point to the conclusion that the Mg II line blend 
is a consequence of the accretion process. The lack of correlation with the accretion rate 
maybe the result of time delay: if the Mg II line is created in the funnel flow one would 
expect a delay of ~ 5 days between a spectral manifestation in Mg II and the landing on 
material on the star. An analogous delay is expected for the wind component of the line. 

On the other hand GdCF97 suggest that magnetospheric flaring may give origin to 
increased Mg II emission. Such flaring is present in the model by Shu et al. (1994). If this 
is the case, it is unclear whether or not any time delay would correlate the line flux with 
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the accretion rate. The flaring of the magnetosphere maybe be causally unrelated to the 
process of accretion. If this is the case we may be seeing a Mg II line corresponding to more 
than one source: a wind component, a funnel component and a magnetospheric component. 

As the model used here to find the accretion rate is sensitive only to regions below the 
accretion shock, the lack of correlation of the Mg II line with the accretion rate is telling us 
that the dominant component of the line is not created in the region below the accretion 
shock. 



5.4. The Accretion rate and the spot size 

Figure 6a shows that there is a correlation between the spot size, /, and the accretion 
rate, M cos9. The slope of the points in Figure 6a is such that Mcos^ a J°- 5±0 - 05 . 

The correlation between Mcos9 and / implies (via Equation 1) a correlation between 
logjF and /, as shown in Figure 6b. Here T a /~ ' 5 . Assuming a constant accretion rate of 
Mcos9 = 1.6 1O~ 8 M yr _1 , we obtain the solid line shown in Figure 6b, which shows that 
smaller spots have under-average Mcos# and larger spots have over-average M cos9. This 
merely reflects the fact that the observed (projected) accretion rate is not constant. 

As Equation 2 shows, the observed (in the detector) ultraviolet flux, in ergs/sec cm 2 
is such that Fjjy ot, Mcos0. This implies, from Figure 6a, that Fjjy a f 05 . The flux on 
the surface of the star is proportional to the detected flux divided by the emitting area. 
Therefore, at the stellar surface, Fjjy a / _0 ' 5 and the particle flux and the ultraviolet flux 
are proportional to each other. 

The proportionality between particle and surface UV flux is implicit in CG98 models 
because a large particle flux produces a large UV flux. What needs an explanation is the 
physical origin of only one of the correlations (Fffy a /~ ' 5 , T a /~ ' 5 or Mcos9 a f ' 5 ). 
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Once one is explained, the others follow. 

CG98 observe values of energy flux and spot sizes suggesting that, for increasing 
mass accretion rate, material will fall onto the star over a larger surface area, while the 
energy flux T remains more or less constant. They reach this conclusion by single optical 
observations of 15 CTTSs (Figure 6c). We also find that a bigger mass accretion rate 
implies a bigger spot size, but the energy flux that we calculate does not remain constant. 
While there is a tendency in CG98's data towards larger accretion spots with larger area 
(as implied by the constant J 7 ), this tendency is not strong. Most of their stars lie on the 
locus defined by the BP Tauri observations and therefore, their data could be construed as 
confirming the analysis done here. However, their stars are not a homogeneous population, 
as they encompass a range of radii and masses. 

The correlations from Figure 6 are somewhat sensitive to the estimate of the 
temperature of the underlying star. A higher temperature star will produce smaller inferred 
accretion spots with higher log J 7 . For the 3580 K, the temperature that we have chosen 
here, the ratio between the observed UV continuum flux and the photospheric flux is ~ 30. 
For 4000 K this ratio ~ 10. For the 4000K case, the accretion spots are 0.3 times the size 
of the spots obtained with the 3580K case. The values of log J 7 and M cos 6 are 1.04 times 
and 0.6 times the values of the 3580K case. The shape of the correlations from Figure 6, 
including slopes and overall values, does not change, within the errors in the calculated 
quantities. In other words, even though a given spot will change size when the assumed 
temperature of the underlying star changes, the overall relationship between spot size, 
energy flux, and accretion rate does not change. 
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5.5. Simulating the Correlation between Accretion Rate and Spot Size 

A possible explanation for the correlation from Figure 6a lies in the geometry of the 
accretion region. As the spot size, /, and the accretion rate, Mcos^, are both proportional 
to the geometric factor, cos 9, one would expect a correlation between them. When the 
projected size of the accretion spot is large (large / due to a small inclination angle) 
we should see large projected accretion rates. The larger range of the spot sizes (~ 1.6 
dex) compared to the accretion rates (~ 1.2 dex) could be explained by different intrinsic 
distributions. 

To decide whether or not the correlation of Figure 6a is due to geometrical effects 
we have perfomed a simulation. We assume that the situation is one in which we are 
looking at an opaque disk from an inclination angle i. The spot has the colatitude angle 
9 co i. As all the spectra in our data set show ultraviolet excess, one always see the spot. 
This gives the constraint that % < n/2 — 6. Let i be the unit vector pointing from the 
center of the star to the observer. Then i = (sin(i), 0, cos(i)). Let n be the unit vector 
perpendicular to the spot area. Then n = (sin(# co ;) cos(atf), sm{9 co i) sin(cjt), cos(6* co z), where 
u = 27r/8 rad per day, the angular frequency of the star, assumed to have a period of 
8 days. Therefore t maps the phase in which we are observing the star. The cosine of 
the angle between the spot and the observer (which we have called here 6) is given by 
cos9 = i • n = sin(i) sin(# co ;) cos{ut) + cos(i) cos(6* co ;). We then generate uniform random 
points between and 8 for t. This procedure should produce the observed distribution of 
cos 6*. 

We assume that M and f int = f / cos 9 have normal distributions. The simulation 
therefore has 6 parameters: the means and sigmas of the normal distributions, the colatitude 
angle of the spot and the inclination of the system. For a given set of parameters, we 
randomly choose 2000 points for M, f int and cos 6* and multiply M by cos 9 and fi nt by the 
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same cos#. The values of f int cos# and Mcos# are constained to be larger than the minima 
of their distributions. In other words, we are assuming that the accretion rate and the 
intrinsic size are uncorrelated and that any correlation comes from the cos 9 that multiplies 
both quantities. 

We do not find any combination of parameters that will satisfy all the constraints 
of the problem: Figures 4a, 4b and 6a. All the simulations that produce a correlation 
between M cos 9 and / give a slope close to 1 for the log-log plot of the two quantities. 
The width of the gaussians determine the scatter of this correlation. This means that the 
basic hypothesis of the model (that M and f int = f / cos 9 are uncorrelated gaussians) is not 
verified by the simulation. 

Dendy, Helander & Tagger (1988) have considered the possibility that the intrinsic 
distribution of accretion rates should be a power law. They arrive to this conclusion based 
on models of self-organized criticality. The range of accretion rates found here is not large 
enough to test this hypothesis. 

While it is clear that a full understanding of Figure 6a should include the effects of 
changing inclinations of the spot as the star rotates, the fundamental explanation is not 
geometry. We have to look for other mechanisms. 

It is possible that some of the inferred changes in energy flux are due to "limb 
brightening" . Assuming that spot size does not change for short cycles, a small spot size 
would imply a high spot inclination, and therefore we would observe higher regions of 
the spot (farther away from the center of the star), which the models predict have higher 
temperatures (Lamzin 1998). Therefore, these higher regions will look like regions of high 
energy flux. One would therefore expect a correlation between log J 7 and /, as observed 
here. If "limb brightening" were indeed an important effect, the slope of the spectra would 
not measure the value of T but the temperature of that layer, which would change for 
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different orientations. The CG98 models assume that the spot is observed from directly 
from above, and therefore it is difficult to account for this effect. Resolving this issue will 
have to wait for more detailed models. 

5.6. Theoretical Expectations of the Relationship between Spot Sizes and 

Accretion Rates 

The total size of the accretion region is determined by geometry of the stellar magnetic 
field. An accretion spot of variable size would therefore imply that the field itself is changing. 
So far, the most sophisticated idea of the accretion process is the so-called 'X-wind' model 
(Shu et al. 1994). This model has been developed for a magnetic dipole aligned with the 
disk, even though most authors agree that such geometry is an oversimplification of the 
more complicated real geometry of the field. In the X-wind model, the fractional coverage 
of the accretion spot with respect to the area of the star (the quantity we have called here 
/) is given by / = y/l - RJR X - V 1 - 1.5 *R*/R X (Ostriker & Shu 1995), where R x is 
the radius of the X-point', that coincides to a few percent with the corotation radius of the 
star and with the truncation radius Rt- Furthermore, R x a M~ 2 / 7 . For R*/Rx ~ 0.2 as 
assumed in this paper, / ~ 6%. Let's assume the accretion rate increases by a factor of 2, 
from the value that produces by R*/Rx ~ 0.2. R x decrease bys a factor of 0.8 and then 
R*/Rx ~ 0.24. Then / ~ 7%. According to Figure 6a, / should increase by a factor of 4 
when the accretion rate increases by a factor of 2. The increase in spot size with accretion 
rate is in the right direction but is not fast enough to explain the correlation of Figure 6a. 

The model by Mahdavi & Kenyon (1998) uses a tilted dipole and consider the accretion 
as occurring onto a thin ring on the surface of the star. For R*/R T ~ 0.2 — 0.3 and a dipole 
inclined by 60 degrees with respect to the plane of the accretion disk, this model predicts 
/ ~ 2 - 4%, and / a M 2 ' 1 which is again a very slow increase in the right direction. 
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It is also clear that the problem remains if we include explicitly the varying truncation 
radius in the correlation from Figure 6a: / a (Mcos#) 2 (l — R^/Rj) 2 . The increase in the 
accretion rate dominates over the slow decrease in the factor £ = (1 — R*/R T ). 

The root of the problem is probably in the assumption of a magnetic dipole as the 
magnetic field of the star. Even though this is a good assumption at large distances from 
the star it is clear that it is an oversimplification if one wishes to explain the geometry of 
the accretion spots, which depends on the field close to the star. 

In Ostriker & Shu (1995) and Mahdavi & Kenyon (1998), the predicted accretion spots 
are on the upper tail of our distribution. This has also been observed by CG98, which 
conclude that the accretion process has to be inhomogeneous. If not all the available area 
is covered at the same time, the observed accretion spot will be a fraction of the available 
area. 

Independently of these considerations one could argue that, according to Figure 6a, 
matter loading on the field lines is not very efficient. We see that the accretion rate increases 
only with the linear dimensions (e.g. the radius or the circumference) of the accretion spot: 
M ~ /°- 5 . So a bigger area produces a more diluted accretion, and therefore, a weaker 
surface flux Fffy. Another way of saying this is that the spot size is very sensitive to 
changes in accretion rate. 

The models that we use here are not complete: they assume a shock parallel to the 
surface of the star, where the real shock may be oblique. They also assume free-fall velocity, 
which, as Ostriker & Shu (1995) show, is not the most appropiate hypothesis. The inclusion 
of these two effects would increase the observed spot size, but it is unlikely that they will 
change the slope of Figure 6a. To rule out this possibility, more complete simulations are 
needed. 
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In summary, the IUE data on BP Tauri puts a strong constraint in any model that 
seeks to explain the size of accretion spots. Accurate modeling should consider the geometry 
of the observations (i.e. inclination angle, spot latitude), a more realistic geometry for the 
magnetic field, oblique shocks, the real velocity of the material falling onto the star, and 
possibly the effects of limb brightening. There may also be an unmodelled efficiency factor 
that determines the way in which flux tubes are actually filled with material, in such a way 
that the accretion rate scales as the radius of the tube. 

6. Conclusions 

We have analyzed all the observations of BP Tauri taken by the International 
Ultraviolet Explorer in the low resolution (AA ~ 6A), long wavelength range (from 
A = 1850A to A = 3350A). We found 61 useful spectra. 

Compared with spectra of standard stars, all the BP Tauri spectra show strong 
ultraviolet excess and color temperatures around 8000 K. Using 43 spectra, we observe 
variability in the ultraviolet continuum of Am cont . ~ 1 magnitude and variability in the 
Mg II line flux of AmM g n ~ 0.8 magnitudes. We observe no correlation between continuum 
flux and Mg II line flux. This result resolves a controversy in the literature regarding the 
relationship of the continuum flux and Mg II line flux in BP Tauri. We conclude that the 
bulk of the Mg II in BP Tauri is not emitted below the accretion shock. Therefore, it may 
be emitted by the accretion column, the stellar wind and/or magnetospheric flaring. We 
do not observe any correlation between color temperature and the average value of the 
continuum. 

Within the context of the magnetospheric accretion paradigm, we have used models 
by CG98 to obtain the energy flux into the star due to the accretion process, and the spot 



-27- 



size for each of the spectra of BP Tauri. We assume a constant inner truncation radius for 
the disk of 5 R*. With 45 spectra, we find an average energy flux of 5 10 11 ergs cm" 2 sec -1 
and an average projected spot size of 4.4 10~ 3 times the area of one hemisphere of the 
star. The latter result coincides with previous observational determinations of spot sizes. 
From the energy flux and the spot size we obtain the accretion rate times cos 9, where 9 
is the average inclination angle of the region that produces the emission. The average of 
this projected accretion rate is 1.6 10" 8 M Q yr" 1 ± 30%, consistent with values found by 
other researchers. This values are weakly sensitive to the value of the underlying star used 
obtain the continuum excess. No correlation is observed between the flux in Mg II line and 
the energy flux, the accretion rate or spot size, affirming the conclusion from the previous 
paragraph. 

The nature of the data and the models is such that we cannot say anything about 
intrinsic variations in the accretion rate or the real (i.e., deprojected) spot size, but we 
conclude that the dispersion in their ratio is M/f ~ 50%, consistent with the errors. 

We observe a correlation between the ultraviolet flux and the accretion rate, similar 
to correlations found by other researchers between the U-band flux and the accretion rate: 
log(L acc /L ) = (1.1 ± 0.1) * log(L ° / L Q ) — (0.5 ± 0.05). This correlation is implicit in 
the models by CG98. Also, we observe a correlation between the accretion rate Mcos# and 
the projected spot size /: Mcos# a y°- 5±0 - 05 . Together, the two correlations imply that the 
UV flux in the surface of the star is proportional to the particle flux. 

These results are significant, in spite of the large errors involved in the determinations 
of each quantity. These correlations stand in contrast to results obtained by Calvet & 
Gullbring (1998). Using individual optical observations of 15 CTTSs, they concluded that 
the accretion rate increases linearly with spot size. A resolution of this discrepancy will 
have to wait until we have applied the magnetospheric accretion models to other stars. 
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Using a simulation, we conclude that the geometry of the system determined by the 
observer, the rotation axis and the inclination of the star can not explain the correlation 
between accretion rates and spot sizes. One must look at other more fundamental 
explanations. 

Models of magnetospheric accretion by Ostriker & Shu (1995) and Mahdavi & Kenyon 
(1998) do predict that the spot size should increase with the accretion rate , but in both 
models, the predicted increase is much slower than what we found here (/ a M 2 ). 

The correlation between accretion rate and spot size could be due to an unmodelled 
efficiency factor that determines how the spot is filled. Conversely, it indicates that the 
spot size is very sensitive to accretion rate. Future modeling should take into account the 
possibility of a non-dipole field for the star, oblique shocks, limb brightening and the effects 
of the geometry determined by the observer. The correlation between accretion rate and 
spot size puts a strong constraint in these models. 

We would like to thank Nuria Calvet and Erik Gullbring for providing us with their 
models. Dr. Calvet in particular spent quite a bit of time helping us to understand 
the nature and strengths of their models. We are grateful to F. Shu and J. Bouvier 
for suggestions that led to the improvement of the text. Thanks are also due to Debi 
Howell-Ardila who edited the manuscript for language. We acknowledge the support of the 
NASA ADP program through grant # NAG5-3471. 
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Figure 1: (a) Average of all the available spectra (solid line) and the normalized 
variance (dotted line). The feature near 3050A is a reseaux mark, (b) Two sample spectra 
(lwpl8976, solid line, lwp09256, dash-dot line) uncorrected by reddening. Notice the strong 
variation in slope and flux level. The position of some spectral features is indicated (Brown 
et al. 1984). 

Figure 2: Mg II line flux vs. UV continuum flux. The continuum flux has been 
calculated from a box 50A wide centered at 2900A. No correlation between the two 
quantities is observed, which means that the Mgll line is not sensitive to the characteristics 
of the accretion spot. 

Figure 3: Average of the logarithm of the flux versus the color temperature. Both 
have been calculated from the continuum fitting procedure (Section 4.2). No correlation is 
observed between the two quantities. 

Figure 4: (a) Histogram for the calculated values of /. / = 4.4 1CT 3 with a dispersion 
of a — 5 10~ 3 . The errors are ~ 10%. The bar indicates the error at / = 0.0095. (b) 
Histogram for the calculated values of log^ 7 , the logarithm of the energy flux into the 
star; logjF = 11.7 and the dispersion is a — 0.2. The error is ~ 1%, as indicated by 
the bar centered in 11.4. With log J 7 and / we calculate Mcos#, the accretion rate 
uncorrected for projection effects. As shown in (c), Mcos6> = 1.6 1O~ 8 M yr" 1 and the 
dispersion is o ~ 0.8 10^ 8 M Q /yr. The errors are ~ 30%, as indicated by the bar for 
Mcos9= 1.3 10- 8 M o yr" 1 . 

Figure 5: The projected accretion rate vs. the UV continuum flux. The UV continuum 
flux is taken from Table 1. This relation is similar to previously found relationships between 
the accretion luminosity and the U-band flux. 

Figure 6: (a) The logarithm of the accretion rate vs. log(/). The boxes indicate the 
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extent of the errors. The correlation is such that M cos(9) a f°- 5±0 - 05 . (b) The correlation 
between logjF and log(/) implied by (a) via Equation 1: T a /~ 5 . The solid line is what 
one would obtain from Equation 1, assuming a constant Mcos^ = 1.6 1CT 8 M yr _1 . (c) 
Same as (a) without the error boxes. The names indicate the positions of stars analyzed by 
CG98 using single optical observations (i.e. AA means AA Tauri, DK means DK Tauri and 
so on for all stars except GM Aur and UY Aur). In their data, some tendency can be seen 
to larger accretion spots with larger accretion rates. 



Table 1. BP Tauri versus standards 



Name 


IUE Name 


Color Temperature a 


Flux Mg II b 


UV Flux c 


BP Tauri 




7000 to 9000 


~ io- 12 


~ io- 12 




lwp09256 


9230 ± 100 




4.15 10~ 12 




lwp 18976 


8380 ± 80 




2.32 10~ 12 


HBC 399 


lwp07540 


3200 ± 500 


8.9 10" 14 


5 10- 14 


61 Cyg. B d 


lwrl2742 


2500 ± 500 


1.2 IO" 15 


1.5 10~ 15 



Note. - - The spectra of BP Tauri and HBC 399 have been corrected for a reddening of 

Ay = 0.51 

a In Kelvins. The standard spectra degrade rapidly to longer wavelengths and so the 
precision of the color temperature is small. 
b In ergs/s cm 2 . 

c Continuum flux in 50A around 2900A. In ergs/s cm 2 . 

d All the quantities in this row have been calculated with the small aperture exposure. The 
fluxes are observed fluxes at 140 pc. 
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Continuum flux at 2900 A (1Q~ 1S ergs/s cm 2 ) 



Table 2. Summary of observations and model independent results 



Name 



RJD 



EqW Mg II 



Flux Mg II 



UV Flux c 



O 

o 
o 



(N 



> 
CO 

m 

CO 

o 
o 
o 

6 
to 

c3 



lwr04192 


3968.36378 


-28.06 


1.4 


2.36 


lwr04192 


3968.36378 


-43.84 


0.8 


0.85 


lwr04249 


3975.56167 


-90.88 


3.2 


1.64 


lwr09309 


4558.89251 






2.61 


lwrlll28 


4809.72996 






1.84 


lwrlll29 


4809.75842 


-117.34 


8.2 


1.98 


lwrlll31 


4809.95337 


-81.04 


2.9 


1.33 


lwrl4424 


5259.90632 


-97.29 


2.5 


1.89 


lwrl4425 


5259.95126 






1.46 


lwrl5003 


5340.37147 


-106.56 


3.9 


1.39 


lwr 15006 


5340.79080 






1.48 


lwr 16945 


5615.85545 


-89.66 


3.0 


1.43 


lwrl6962 


5617.98471 


— 


— 


2.02 


lwr 169 72 


5620.98243 


-98.44 


2.9 


0.96 


lwrl6980 


5622.98617 


-49.99 


1.9 


1.12 


lwp09227 


6705.42338 


-83.74 


3.1 


1.27 


lwp09228 


6705.55286 


-78.76 


2.9 


1.34 


lwp09229 


6705.73196 






1.48 


lwp09230 


6705.81649 






1.73 


lwp09231 


6705.86626 


-101.77 


2.9 


1.55 


lwp09254 


6709.40975 


-36.39 


4.1 


4.17 


lwp09255 


6709.53838 


-44.33 


4.1 


4.43 


lwp09256 


6709.73180 






4.15 



Table 2 — Continued 



Name RJD EqW Mg II a Flux Mg II b UV Flux c 



lwp09257 


6709.81775 


-40.30 


3.3 


3.58 


lwp09258 


6709.85730 


-39.90 


3.9 


3.78 


lwp09283 


6713.75420 


— 


— 


2.42 


lwp09284 


6713.81044 


— 


— 


2.10 


lwp09285 


6713.86425 


— 


— 


1.98 


lwp09336 


6720.42407 


-72.68 


2.8 


1.59 


lwp09337 


6720.68121 


— 


— 


1.45 


lwp09375 


6724.43171 


-64.49 


3.9 


2.82 


lwp09376 


6724.68419 


— 


— 


4.15 


lwp09416 


6730.43992 


-87.49 


2.4 


1.08 


lwp09418 


6730.68940 


-201.9 


4.4 


1.11 


lwpl4276 


7454.00667 


— 


— 


1.94 


lwpl4791 


7535.54378 


— 


— 


2.10 


lwpl6660 


7825.82212 


-55.30 


2.2 


1.51 


lwpl7185 


7911.53778 






2.75 


lwp 18976 


8174.79117 






2.32 


lwp22194 


8627.20920 


-73.96 


3.9 


1.99 


lwp22195 


8627.44427 


-79.56 


4.1 


2.19 


lwp22211 


8629.18550 


-77.77 


5.2 


2.73 


lwp22212 


8629.24487 


-98.36 


5.7 


2.28 


lwp22213 


8629.44565 


-104.68 


5.1 


2.60 


lwp22226 


8631.17911 


-90.45 


3.4 


1.55 


lwp22227 


8631.30942 


-70.53 


3.4 


2.00 



Table 2 — Continued 



Name 


RJD 


EqW Mg II a 

1 o 


Flux Mg II b 


UV Flux c 


lwp22228 


8631.44646 


-77.83 


3.5 


1.92 


lwp22238 


8633.16363 


-114.15 


5.0 


1.95 


lwp22239 


8633.28937 


-92.46 


4.7 


2.16 


lwp22240 


8633.44058 


-103.97 


4.5 


1.87 


lwp22243 


8635.17840 


-79.89 


3.7 


1.45 


lwp22244 


8635.30686 


-134.56 


4.5 


1.40 


lwp22245 


8635.44097 


-189.39 


4.4 


1.25 


lwp22248 


8637.17266 


-119.27 


3.2 


1.37 


lwp22249 


8637.44297 


-77.94 


4.9 


2.91 


lwp22261 


8639.18842 


-70.00 


4.0 


2.38 


lwp22262 


8639.31788 


-71.68 


3.9 


2.59 


lwp22263 


8639.44547 


-81.71 


4.2 


2.41 


lwp22276 


8641.17006 


-58.75 


3.6 


2.11 


lwp22277 


8641.30575 






1.40 


lwp22278 


8641.44384 


-78.94 


3.9 


1.84 



Note. - - The empty positions indicate spectra in which the "quality" flags for pixels in 
the Mg II lines are different from zero. 

Equivalent Width of the Mg II blend. Units: A. The error is ±10A 

b Net flux in the Mg II blend. Units: 10" 12 ergs/s cm 2 . Error: 0.05 10 12 ergs/s cm 2 

c Flux in a 50A box around 2900A. Units: 10" " 12 ergs/s cm 2 . Error: 0.02 10 12 ergs/s cm 2 
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Color temperature of the continuum 
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o 

O ' Table 3. Modeling the UV continuum 



(N 



> 
CO 

m 
m 
o 
o 
o 

a: 

6 
to 

c3 



13 



Name 




Slope b 


Average 




log.F e 


M cos l 


l r\ a i no 

lwr04192 


7onn / a r\r\\ 

7300(400) 


O *7 A 1 l '\ — 4- 

2.74 10 


1 o r - no 

—13.502 


O 1 / 1 1 \ in — 

8.1(1.1) 10 J 


1 i a { n n \ 

11.4(0.2) 


1 O ( 1 n \ 

1.8(1.0) 


l n /i i no 

lwr04192 














1 n a o /1 n 

lwr04249 














lwr09309 


8700(500) 


1.03 10 4 


-13.448 


2.5(0.3) 10 6 


11.8(0.1) 


1.6(0.5) 


lwrlll28 














lwrlll29 














lwrlll31 














lwr 14424 














lwrl4425 


7800(900) 


1.38 10~ 4 


-13.752 


1.6(0.3) 10~ 3 


11.7(0.3) 


0.8(0.6) 


lwrl5003 














lwrl5006 


6900(300) 


3.18 10~ 4 


-13.738 


7.9(0.8) 10- 3 


11.2(0.2) 


1.2(0.6) 


lwrl6945 


7500(800) 


1.66 10~ 5 


-13.787 


1.8(0.3) 10~ 3 


11.6(0.3) 


0.8(0.6) 


lwrl6962 


9700(900) 


0.00 


-13.589 


9.4(2) 10~ 4 


12.0(0.2) 


1.0(0.5) 


lwr 169 72 














lwrl6980 














lwp09227 


8000(500) 


9.09 10~ 5 


-13.835 


9.8(1) HT 4 


11.8(0.2) 


0.6(0.2) 



Table 3 — Continued 



Name 




Slope b 


Average 


/ d 


o 


M cos 6> f 


lwp09228 


8100(500) 


9.20 10~ 5 


-13.790 


1.1(0.1) 10~ 3 


11.8(0.1) 


0.7(0.2) 


lwp09229 


7900(200) 


1.87 10~ 4 


-13.638 


2.5(0.2) 10~ 3 


11.6(0.1) 


1.1(0.2) 


lwp09230 


9700(900) 


-2.58 10~ 5 


-13.707 


6.2(1) 10~ 4 


12.1(0.2) 


0.8(0.3) 


lwp09231 


- 


- 


— 


— 


— 




lwp09254 


8301(200) 


1.80 10~ 4 


-13.157 


7.2(0.5) 10- 3 


11.7(0.1) 


3.4(0.5) 


lwp09255 


9100(200) 


9.28 10~ 5 


-13.160 


4.1(0.2) 10- 3 


11.9(0.1) 


3.0(0.4) 


lwp09256 


9230(100) 


1.18 10~ 4 


-13.170 


3.5(0.2) 10" 3 


11.9(0.1) 


3.0(0.2) 


lwp09257 


7600(200) 


2.61 10~ 4 


-13.305 


1.0(0.08) 10~ 2 


11.4(0.1) 


2.8(0.7) 


lwp09258 


7800(200) 


2.41 10~ 4 


-13.200 


1.1(0.07) lO" 2 


11.5(0.1) 


3.4(0.7) 


lwp09283 


8500(400) 


1.10 10~ 4 


-13.518 


2.2(0.2) 10" 3 


11.8(0.1) 


1.4(0.3) 


lwp09284 


8900(200) 


9.89 10~ 5 


-13.472 


2.1(0.1) 10~ 3 


11.8(0.1) 


1.5(0.2) 


lwp09285 














lwp09336 


7400(400) 


2.16 10" 4 


-13.708 


3.2(0.3) 10~ 3 


11.5(0.2) 


1.0(0.4) 


lwp09337 


7200(300) 


2.51 10~ 4 


-13.716 


4.3(0.3) lO" 3 


11.4(0.2) 


1.1(0.4) 


lwp09375 


7400(200) 


2.66 10~ 4 


-13.365 


1.0(0.06) 10~ 2 


11.4(0.1) 


2.5(0.7) 


lwp09376 


7200(100) 


3.37 10~ 4 


-13.221 


2.6(0.1) 10~ 2 


11.2(0.1) 


4.2(0.7) 



Table 3 — Continued 



Name 




Slope b 


Average 


r 






\ogF c 


Mcos6» 


lwp09416 


- 


- 


— 


— 






— 




lwp09418 


— 


- 


- 


- 






- 


— 


lwpl4276 


8900(200) 


8.14 10~ 5 


-13.490 


1.9(0.1) 


10" 


-3 


11.7(0.1) 


1.4(0.2) 


lwpl4791 


8700(100) 


1.33 10~ 4 


-13.465 


2.5(0.1) 


10" 


-3 


11.8(0.1) 


1.6(0.1) 


lwpl6660 


- 


- 


— 


— 






— 




lwpl7185 


9600(200) 


4.16 10~ 5 


-13.373 


1.8(0.1) 


10" 


-3 


12.0(0.1) 


1.7(0.1) 


lwpl8976 


8400(200) 


1.58 10~ 4 


-13.437 


3.2(0.1) 


10" 


-3 


11.7(0.1) 


1.7(0.2) 


lwp22194 


7300(300) 


2.62 10~ 4 


-13.539 


6.5(0.5) 


10" 


-3 


11.4(0.1) 


1.6(0.5) 


lwp22195 


9000(400) 


6.24 10~ 5 


-13.528 


1.6(0.1) 


10" 


-3 


11.8(0.1) 


1.3(0.2) 


lwp22211 


8600(400) 


1.06 10~ 4 


-13.457 


2.5(0.2) 


10" 


-3 


11.8(0.1) 


1.6(0.3) 


lwp22212 


9500(800) 


-1.75 10^ 5 


-13.540 


1.1(0.1) 


10" 


-3 


12.0(0.1) 


1.1(0.4) 


lwp22213 


7800(400) 


1.88 10~ 4 


-13.469 


4.4(0.5) 


10" 


-3 


11.6(0.1) 


1.7(0.6) 


lwp22226 


8600(600) 


2.65 10~ 5 


-13.818 


7.1(1) 


10" 


4 


11.9(0.2) 


1.7(0.2) 


lwp22227 


8000(300) 


1.79 10^ 4 


-13.530 


3.3(0.2) 


10" 


-3 


11.6(0.1) 


0.6(0.3) 


lwp22228 


7800(300) 


2.10 10~ 4 


-13.575 


3.6(0.3) 


10" 


-3 


11.6(0.1) 


1.4(0.3) 


lwp22238 


7300(400) 


2.63 10~ 4 


-13.625 


5.7(0.7) 


10" 


-3 


11.4(0.2) 


1.4(0.7) 



Table 3 — Continued 



Name 




Slope b 


Average 




log.F e 


M cos e { 


lwp22239 


9400(400) 


2.24 10~ 5 


-13.535 


1.2(0.1) 10~ 3 


11.9(0.1) 


1.2(0.2) 


lwp22240 


8600(400) 


1.03 10~ 4 


-13.591 


1.8(0.1) 10~ 3 


11.8(0.1) 


1.1(0.3) 


lwp22243 


8600(500) 


5.73 10~ 5 


-13.693 


1.1(0.1) 10~ 3 


11.8(0.1) 


0.9(0.3) 


lwp22244 


- 


- 


- 


- 


- 


- 


lwp22245 


8000(500) 


8.39 10~ 5 


-13.803 


1.2(0.1) 10~ 3 


11.7(0.2) 


0.7(0.3) 


lwp22248 


8300(500) 


1.17 10~ 4 


-13.755 


1.2(0.1) 10~ 3 


11.8(0.2) 


0.8(0.3) 


lwp22249 


7200(200) 


3.24 10~ 4 


-13.329 


1.6(0.1) 10~ 2 


11.3(0.1) 


3.0(0.9) 


lwp22261 


8200(300) 


3.00 10~ 4 


-13.467 


3.4(0.3) 10~ 3 


11.7(0.1) 


1.6(0.3) 


lwp22262 


8900(300) 


8.56 10~ 5 


-13.454 


2.1(0.1) 10~ 3 


11.8(0.1) 


1.5(0.3) 


lwp22263 


8400(300) 


1.38 10~ 4 


-13.486 


2.7(0.2) 10~ 3 


11.7(0.1) 


1.5(1.2) 


lwp22276 


7000(300) 


3.24 10- 4 


-13.557 


1.5(0.01) 10~ 2 


11.1(0.2) 


1.9(0.3) 


lwp22277 














lwp22278 


8700(400) 


8.28 10~ 5 


-13.642 


1.4(0.1) 10~ 3 


11.8(0.1) 


1.0(0.2) 



Note. - Empty positions in the table indicate spectra that are so noisy that a 

determination of the continuum over a sizable wavelength range is impossible. Numbers 

in parenthesis indicate errors 
a Color temperature in K 

a Slope of the logarithm of the fitted flux. 

b Average of the logarithm of the fitted flux. 

c Size of the accretion spot as fraction of the observed area of the star. 

d Logarithm of the energy flux 

c Proiected accretion rate. Units: 10~ 8 M^/vr 
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